Can amphiphile architecture directly control vesicle size? 
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Bilayer membranes self-assembled from simple amphiphiles in solution always have a planar 
ground-state shape. This general principle is a robust consequence of several internal relaxation 
mechanisms of the membrane and prevents the straightforward control of the size of vesicles. Here, 
we show that this principle can be circumvented and that direct size control by molecular design is 
a realistic possibility. Using coarse-grained calculations, we design tetrablock copolymers that form 
membranes with a preferred curvature, and demonstrate how to make the free energies of these 
structures lower than those of other self-assembled geometries. 



A fundamental process in soft matter science is the 
self-assembly of amphiphilic molecules into a startling va- 
riety of structures, ranging from simple micelles to com- 
plex connected aggregates Self-assembled struc- 
tures not only occur naturally in living cells, but can 
also be designed for vital applications such as drug en- 
capsulation and delivery Q . The question that lies at the 
heart of this field is how the properties of the individual 
amphiphilic molecules control the topology of the aggre- 
gates they form || . One of the major unsolved problems 
is whether a molecule can be designed that can directly 
fix the curvature of a membrane in solution. In addi- 
tion to its fundamental interest, this question is of great 
practical importance, as finding such a molecule would 
allow the straightforward and spontaneous formation of 
vesicles of a well-defined size, thereby yielding precise 
control of drug delivery systems. 

At present, membrane curvature can only be controlled 
by rather complex procedures. Several of these lo HllI re- 
quire two species of amphiphile to be blended |12j |. so 
that the symmetryof the inner and outer leaflets of the 
bilayer is broken [13( and the vesicle has a preferred ra- 
dius. Such methods have the disadvantage that blends 
of amphiphiles can form a wide range of micellar aggre- 
gates which may coexist with the target vesicle struc- 
ture Q . Other methods involve the use of more compli- 
cated vesicle formati on p athways, such as dewetting from 
a templated surface [14J , cooling and warming through a 
cylinder- vesicle shape transition [l5| , electroformation on 
micropatterned glass slides , flow focusing in microflu- 
idics [IT}, and combined extrusion and dialysis fljj ]. 

In this Letter, we investigate an alternative strategy 
for controlling membrane curvature. We break the mem- 
brane symmetry by the use of ABCA' tetrablock copoly- 
mers [19h21| ]. The outer A and A' blocks of the poly- 
mer are formed of the same hydrophilic material, and 
the B and C blocks are hydrophobic and have a repul- 
sive interaction with each other. These molecules form 
asymmetric monolayers in solution [l9j (see Fig. [TJi), in 
contrast to the bilayers formed by diblocks (Fig.QJi). We 
use tetrablocks rather than ABC triblocks since, in this 
latter case, the A and C blocks would have to be hy- 



FIG. 1: (a) Bilayer vesicle formed from diblock copolymers. 
Hydrophilic blocks are light gray and hydrophobic blocks are 
dark gray, (b) Monolayer vesicle formed from tetrablock 
copolymers. There are now two hydrophobic blocks, colored 
black and dark gray respectively. The black block is more 
strongly hydrophobic, and segregates to the inner half of the 
monolayer in order to minimize the area of its interfaces with 
other species. 



drophilic and have a strong mutual repulsion for asym- 
metric monolayers to form. This combination is hard 
to achieve, both because of the difficulty of finding hy- 
drophilic compounds that repel strongly and the fact that 
the hydrophilic layers are diluted by solvent, weakening 
any interaction between them. 

Tetrablock vesicles have indeed been formed in prelim- 
inary experimental investigations by Brannan and Bates 
19], who also achieved a degree of size control of the 
aggregates. Such vesicles have also been found in very 
recent Monte Carlo simulations 21 1, although the ques- 
tion of size control was not studied. Here, we present 
a concrete theoretical demonstration of the basic prin- 
ciple that tetrablock copolymers can form bilayers of a 
preferred curvature, and also show how these molecules 
should be designed to control the vesicle radius and avoid 
the formation of unwanted micellar structures [20j ]. 

To capture the basic physics of the system as clearly 
as possible, we focus on a simple model of ABCA' tetra- 
blocks blended with A homopolymer 'solvent'. To begin, 
we consider copolymers where all four segments contain 
the same number of monomers N/4. For simplicity, the A 
homopolymer molecules were also taken to contain N/4 
monomers. The strength of the interactions between the 
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various species are set by Flory x parameters. We note 
that once two x parameters have been chosen, the third 
must be calculated from a relation involving the polar- 
izabilities of the various species [23 |. To calculate the 
density profiles and free energies of the self-assembled 
structures, we used a simple coarse-grained mean- field 
theory (self-consistent field theory, or SCFT) in which 
the individual polymer molecules are modeled by ran- 
dom walks and their interactions by a contact potential 
H Hit The great advantage of SCFT for our current 
investigation is that it is simpler and faster than, for ex- 
ample, Monte Carlo methods, allowing the free energy 
to be computed for a large range of vesicle sizes and 
the question of optimum bilayer curvature to be directly 
addressed. Furthermore, for the high molecular weight 
synthetic polymers considered here, SCFT can approach 
Monte Carlo simulations in accuracy [25[ and can yield 
considerable qualitative insights even into aqueous so- 
lutions of these molecules [9(. The diffusion equations 
describing the random-walk polymers were solved by a 
standard finite-difference method and the SCFT equa- 
tions by an iterative scheme (26| . 

Since we focus on curved membranes, and in particular 
on spherical vesicles, we perform many of our calculations 
assuming spherical symmetry in a spherical calculation 
box, although other geometries will also be considered. 
As these calculations involve only a single aggregate, we 
must make the connection between the free energy of 
the subsystem of volume V containing the vesicle and 
that of the system as a whole. To do this, we calculate 
the free-energy density F of a box containing a single 
spherically-symmetric vesicle in the canonical ensemble. 
We measure all energies in units of k^T, and make F 
dimensionless by multiplication by the volume v of a sin- 
gle copolymer molecule. We then vary the volume of 
the simulation box at constant overall copolymer volume 
fraction [27j |. In a very simple way, this mimics a larger 
system (of fixed total volume and fixed copolymer vol- 
ume fraction) searching for its free energy minimum by 
varying the number of aggregates and hence the volume 
occupied by each. If the calculated free energy density 
has a minimum as a function of the volume of the calcula- 
tion box, this means that the aggregates have a preferred 
size, and, in our case, that the vesicle membranes have 
an optimum curvature. 

Although this minimum corresponds to the absolute 
free energy minimum of a solution of spherical vesicles, 
the shape of the curve F(V) is not directly linked to the 
size distribution of the aggregates in a real system, as a 
point on the curve corresponds a monodisperse solution 
of vesicles of a given size. Fortunately, these results can 
be interrogated further to calculate the polydispersity A 
of the vesicles, which arises from fluctuations around the 
free energy minimum. In theories of micellization [28j j . 
this is related to the curvature of the free energy f p of 
an aggregate of p molecules by 1/A 2 = d 2 f p /dp 2 . To ex- 
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FIG. 2: Sample volume fraction profiles for (a) ABCA' vesi- 
cles, (b) vesicles in the unfavorable A'CBA orientation, and 
(c) micelles. A-blocks: thin full lines, B-blocks: thick full 
lines, C-blocks: thick dashed lines, A'- blocks: thin dashed 
lines. 



tract this quantity from our curves of F{V), we write the 
(dimensionless) free energy density of our monodisperse 
solution of aggregates in terms of f p as 

F = (<f>-v m /V) \n[{<j)-v m /V)/e] + {<t>-v m /V)h+vf p /V 

(1) 

where <fr is the volume fraction of copolymer molecules, 
1/V is the number density of aggregates, v m — pv is 
the volume of an aggregate, /i is the energy of an iso- 
lated copolymer in solution, and f p is the energy of an 
aggregate of p copolymers. The first term arises from 
the entropy of the free copolymers in solution. Now, we 
note that a single SCFT calculation finds the local free 
energy minimum F(V) for a vesicle in a box of volume 
V. In the process, it determines the optimum number 
of molecules in the vesicle for this box size and there- 
fore corresponds to minimizing F with respect to p at a 
given 1/V. Varying the box size then yields the curve 
F(V), from which we can read off the second derivative 
d 2 F/dV 2 . Remembering that this derivative is evaluated 
along the line where dF jdp\y = 0, we find, by differen- 
tiation and rearrangement of Eq. [TJ that 
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(2) 



allowing us to calculate A. 

We now demonstrate the existence of the target vesicle 
structure as a solution to SCFT. In Fig. [2p, we plot cuts 
through the density profile of an ABCA' vesicle, with 
the overall copolymer volume fraction fixed to 0.05. The 
various x parameters must be large enough for the am- 
phiphile to aggregate, and so we set xab = 50/A and 
Xac — 30/iV, where N is the total number of monomers 
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FIG. 3: Minimum of the free energy density as a function of 
vesicle radius R. The free energy Fh of a homogeneous system 
of the same composition is subtracted and the free energy den- 
sity is made dimensionless by multiplication by the volume v 
of a single molecule. The inset shows the corresponding curve 
for vesicles formed of tetrablock polymers in the unfavorable 
A'CBA orientation. 

in the copolymer. So that the B and C species demix, 
we choose the larger of the two possible values [22[ for 
Xbc, which we set to 157. 5/N. We measure all leng ths 
in units of a_/V 1//2 , where a is the segment length [24|. 

The ABCA' structure sketched in Fig. [T]is clearly re- 
produced in our calculations. The strongly hydrophobic 
B-blocks lie in the inner half of the membrane, so that 
the more energetically unfavorable AB interface has a 
smaller area than the A'C interface. Figure [3^ shows 
that the local density is highest in the B layer, and lower 
in the more diffuse hydrophilic corona. We note also that 
the density profiles of both hydrophilic blocks have small 
peaks in the center of the vesicle wall. This is because 
a fraction of the copolymers fold back on themselves to 
protect the strongly incompatible B and C species from 
each other. Surprisingly, we also find solutions to SCFT, 
shown in Fig. [2p, where the vesicle is 'inside-out', with 
the B species on the outside of the vesicle core, increasing 
the size of the AB interface. 

In Fig. we plot the density profile of a micelle, 
formed in a smaller calculation box. Here, the micelle 
core is formed from the strongly hydrophobic B blocks, 
while the other species mix in the corona. This structure 
is most likely formed as in ABA triblock micelles [29j j . 
with the copolymers forming a hairpin shape. 

We now focus on the free energies of the structures 
plotted in Fig. [21 and in particular on their size depen- 
dence. In Fig. [31 we plot the free energy density Fasa 
function of the ABCA' vesicle radius R (defined as the 
radius on the outside of the vesicle at which the copoly- 
mer and solvent densities are equal), which we vary by 
changing the calculation box size as detailed above. For 
concreteness, we fix the zero of our free energy scale 
to correspond to a homogeneously-mixed system of the 



same composition. This choice does not affect our cal- 
culations of polydispersity, which require only the cur- 
vature d 2 F/dV 2 . The curve shows a clear minimum as 
a function of R, demonstrating that the vesicles have a 
preferred size. This in sharp contrast to the the mono- 
tonic decrease of the free energy density of the A'CBA 
vesicle, plotted in the inset. 

An insight into the mechanism behind the existence of 
an optimum vesicle size can be gained by comparison of 
the two curves. At smaller radii, the free energy den- 
sity of both vesicles decreases with increasing size, as the 
chains relax and the copolymers are no longer strongly 
compressed on the inner face of the membrane. The 
ABCA' vesicle always has a lower free energy than the 
A'CBA vesicle, as the unfavorable AB interface is on the 
inside and is therefore smaller. As the vesicle radius in- 
creases, both membranes become flatter and the relative 
advantage of the ABCA' vesicle decreases. The two lines 
then start to approach each other, with the ABCA' curve 
now rising and the A'CBA curve continuing to fall. The 
monotonic form of the A'CBA curve is also clear evidence 
that the ABCA' minimum is not a finite size effect. 

To calculate the relative polydispersity A/p of the vesi- 
cles, we plot the main free energy curve of Fig. [3] as a 
function of subsystem volume V , allowing us to calcu- 
late the curvature d 2 F/dV 2 at the minimum. We cal- 
culate the aggregate volume v m by integrating over the 
vesicle density profile and subtracting the local volume 
fraction at the edge of the system, where it has attained 
a stable bulk value. To proceed further with this calcu- 
lation, we estimate the volume v of a single copolymer 
molecule. By recalling that all volumes are measured in 
units of a 3 N 3 / 2 , and defining the segment volume such 
that v = a 3 N, we can show that A/p is given by the 
product of a term specified uniquely by our SCFT cal- 
culations and l/N 1 / 4 . This shows that the polydisper- 
sity is rather insensitive to the choice of degree of poly- 
merization within the physically reasonable range [20| of 
N ~ 100 — 1000. Even using the smallest of these values, 
N = 100, we find clear size selection, with a relative poly- 
dispersity of A/p w 0.09. Since the vesicle is relatively 
flat, we can assume that its surface area is proportional 
to p and hence that its radius is proportional to p 1 / 2 . 
This yields a relative polydispersity in terms of the ra- 
dius of 0.05, and shows that strong size selection can be 
obtained in our simple model system. 

By reducing the simulation box size further, we also 
find a minimum in the free energy density corresponding 
to a spherical micelle (not shown in Fig. [3]). However, 
for the strongly repulsive interactions studied here, this 
structure, in which the incompatible A and C species 
must mix in the corona (see Fig. [2J;) , has a much higher 
free energy density than the ABCA' vesicle, with F « 
—0.034. By performing calculations in a cylindrical box 
at the same copolymer volume fraction, we can also find a 
free-energy minimum corresponding to an (infinite) cylin- 
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drical micelle, which, at F as —0.037, lies in between 
those of the spherical micelle and the ABCA' vesicle. The 
ABCA' vesicle therefore has the lowest free energy of all 
simple structures for the current parameters. This sit- 
uation changes if the hydrophobicity of the C-blocks is 
reduced so that xac = 20/N and xbc = 133.2/JV. In 
this case, the spherical micelle has the lowest free energy 
density, with F as —0.028, since the A and C species can 
now mix more freely in the corona. The next lowest free 
energy is that of the cylindrical micelle, F as —0.023, 
while the free energy of the ABCA' vesicle is now much 
higher than that of the micelles, with F as —0.0076. Re- 
cent experiments have indeed shown that ABCA' copoly- 
mers can form micelles rather than planar structures at 
room temperature , and this result demonstrates how 
copolymers with sufficiently high might be used to 
avoid this outcome. 

To address the important practical question of how the 
architecture of the copolymer affects the optimum vesi- 
cle radius, we repeat the calculation of the free energy 
minimum shown in Fig. [3] with the size of the B-blocks 
increased by 10%. This leads to a growth of the radius 
from R as SAaN 1 / 2 to R as QAaN 1 ' 2 , to avoid unfa- 
vorable compression of the longer B-blocks on the inner 
half of the vesicle membrane. This opens the possibility 
of rather fine control of the vesicle size by the degree of 
polymerization of the hydrophobic blocks. 

In summary, we have demonstrated that the curvature 
of membranes in solution can be controlled by the archi- 
tecture of the constituent amphiphilic molecules. Specif- 
ically, we use coarse-grained calculations to show that 
copolymers composed of two central hydrophobic blocks 
and two outer hydrophilic blocks form vesicles with a pre- 
ferred radius. To our knowledge, this is the only system 
where the molecular structure of the amphiphiles can be 
shown directly to fix the curvature of a membrane in 
solution. Control of curvature has only been achieved 
before by mixing two types of amphiphile or by using 
a mechanically complex self-assembly method. We have 
also shown how to encourage the formation of vesicles by 
the use of appropriate block chemical nature and how to 
adjust the preferred radius by changing the block length. 

Several practically important extensions to our re- 
search suggest themselves. Our current coarse-grained 
methods can be implemented in a straightforward man- 
ner to explore the chemical parameter space more com- 
prehensively; in particular, the role of the various block 
lengths and interaction strengths in fixing vesicle size 
and polydispersity could be investigated fully. Also, the 
model can be strengthened by introducing further molec- 
ular details, for example using dissipative particle dy- 
namics [13] or Monte Carlo (2]J, with a particular view 
to realistic modeling of the solvent and to optimizing the 
experimental polymer parameters for the production of 
nearly monodisperse vesicles. 
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